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ABSTRACT 

Because of the ccxnplexity of most three- and four-center integrals that 
arise in the variational treatment of molecular structure, approximations have 
been useful in their evaluation. The Margenau approximations are intended to 
have good accuracy when the distance, ft, between the charge distributions involved 
in the integrals is large compared to the spatial size of the distributions. 
Comparisons of values given by these approximations with accurate values of 
some relatively simple three- and four-center integrals confirm this fact. 
Comparisons with the Sklar and Mulliken approximations indicate that Margenau 's 
second approximation is more accurate than those latter approximations in nearly 
1 cases similar to those considered here when ft i*5 at least about five times 
the raidial size of the atomic orbitaG.s upon which the molecular wave functions 
are based. The form of the multipole expansion for the integrals may be used 
to predict euad Interpret the relative accuracies of the approximations in terms 
of the inclusion or neglect of the effects of higher multipoles and cverlaping 
of the charge distributions. A proposed combination of Margenau 's two approxi- 
mations is found to be superior to both in most cases. 


This work was supported in part by the Office of Naval Research and the 
Air Force Office of Scientific Research (at Yale University). 

^A major portion of this work was carried out while the author was an 
Instructor in Physics at Yale University. 

^^The author is a National Research Council Postdoctoral Associate in the 
Laboratory for Space Physics at Goddard Space Flight Center. 



I. HWRODUCTION 


In spite of the large amount of work that has been done on the evaluation 

of the three- euad four-center integrals that occur in the variational treatment 

of molecular structure, meuny of these integ^’als remain a problem because of 

their great complexity. Except in a few cases the exact forms which have been 

found for some of the integrals involve infinite series and other expressions 

or procedures difficult to use.^ Although computers have been successfully 

2*4 

employed in the numerical evaluation of some of the integrals , the results 

often provide only limited insight into the nature of the integrals. 

Several approximations for the three- and four-center molecular integrals 

5“17 

have been proposed and have been used with apparent success. Consideration 

5 

is given here to two approximations proposed by Margenau which are not spec- 
ifically well known, but which have yielded good I'esults for the H-Ha and 

5-7 

intermole cular potential energies. 

Direct determinations of the accuracy of these approximations by competrisons 

of accurate euad approximate values are presented here for a number of relatively 

simple three- and four-center molecular Integrals Involving certain specific 

nuclear configurations. In addition, their accuracies are compared with the 

8 9 

accuracies of the Sklar and Mulliken'^ approximations. The results are capable 
of being predicted euad interpreted in terms of the degree of neglect or over- 
emphasis of the effects of overlap and higher multipoles of the charge distri- 
butions involved. 
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II. TBE THREE- AHD FOUR-CEHTER MOLECULAR INTEGRALS EMPLOYED 
The Integrals for which comparisons oetweeu accura -c o-d aj-proximate veJ.ues 
are given in Sec. VI are instances of the three- and four-center coulomb exchange 
Integrals, 

where ¥ , Y. Y . Y., are normalized real atomic wave functions whose forms 

aor bp cy do 

are designated by a, 0, y, 6 and which are centered on nuclei at positions 

a, b, c, d; rx and rg are general designations of the coordinates of elecurons 

1 and 2} and rxs is the interelectron distance. 

Figure 1 show’ «.he -'axious nuclear centers, distances, emgles, line seg- 

ment% and points involved in the discussions of this and following sections. 

Each integral given by Eq. (l) is equivalent to the electrostatic potentisil 

energy between two distributions of charge, eY„ (ri) Y, o(rx) and eY . Gz) Y, Gz). 

'aor bP CY dy 

The integrals are four-center integrals when the four nucleeur centers a, b, 
c, d are distinct, three-center integrals when two of the centers coincide. 

In the latter case the resulting Integrals are the absolute values of the negative 
nuclear attr8w:tion integreJ.s if the cheurge distribution about the coinciding 
centers (taken here to be c and d) is a delta function of charge e. Hence, 

Eq. (1) may be considered to include the nuclear attraction integrals that are 
likewise considered in Sec. VI. 

The integrals given by Eq. (l) are the only integrals involving more than 
two centers which arise in the variational treatment of molecular structure, 
provided the wave function for the molecular electrons is taken to be a 
linear combination of products of atomic electron wave functions (LCAO method). 
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The Integrals considered here result from Eq. (l) vhen 

Y. o(ri) are each taken to be the ground state hydrogen wave functions for 

unit nuclear charge, (n a ^ e and (tt a ^ e while 

' o o ' 

Y^^(ra) and either taken to be of the same form, or the charge 

distribution e of charge e. Each of the 

resulting Integrals then depends only on the positions of the centers 
a, b, c, d. 

The approximations and accurate- values are compared for two types of 
sequences of positions of the centers. In the first type, the distance, R, 
between 6uad which sure the midpoints of line segments ab and cd, is 
varied while the lengths, and of line segments ab and cd are fixed. 

In the second type of sequence the ratios Rgjj/® aiid ®cd^^ fixed while R 
varies. The sequences include the limits R -♦ 0 and R 

In the first type of sequence and in the limit R -• 0 some three- and 
four-center integrals remain difficult to evaluate, whereas others, like 
those employed in Sec. VI, reduce to slnqpler forms which may be calculated. 

In the limit R -• « all integrals represented by Eq. (l) approach Qx Qa/'* 
this sequence. The quantities Qx and Qg are the total charges of the 
distributions of electrons 1 and 2. In the second type of sequence and in the 
limit R -• 0 all integrals giver by Eq. (IJ become single center integrals, 
the values of which are known or may be computed. For R • in this sequence 
the limiting forms of the integreds are not always easy to evaluate. In this 
limit most, if not all, integrals for which lines ob and cd overlap for a portion 
of their lengths have a form, const, x Qx Qe JJn R/fe, whereas most, if not all, 
of the other integrals reduce to const, x Qx 
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III. A MEANS FOR JUTERERETING THE ACCURACIES 

t' ' 

The bipolar multipole expansion of the integrals given by Eq. (l) may be 
used in the qusJ.itative consideration of the nature and accuracies of the 
Margenau approximations as vel^ as in comparing them with the Sklar and Mulliken 
approximations. The centers for the bipolar expansion are chosen to be ab 
and since these points are the geometric centers of the charge distributions 
involved in the Integrals considered here. The bipolar multipole expansion 
may be obtained by expanding the Integrsds in spherlced hsurmonlcs of the polsur 
angles 6j. and @2 between lines ab and cd and the line aScd, and the azimuthal 
angles ix lines ab and cd around line Sb^. The result as applied to 

Eq. (l) may be written in the fomj, 

"A 


• UmV COSM(i-4) , 


( 2 ) 


|mi 


where the are associated Legendre polynomials, while the symbols > 

auad < denote 'the greater of and 'the lesser of respectively. In the type 
of multipole eiqpanslon represented by Eq. (2) the Is^llclt multipole moments 
aire functions of R in addition to and and the expansion is exact for 
all values of R, R^ aiad R^^. 

A term in the expansion, (2), when the sum over n» has been carried out, is 

th 

the contribution to I from the potential energy between the tx multipole 

th 

of the charge distribution of electron 1 and the is multlpol%of the 
charge distribution of electron 2. Because of the overlap of the charge 
distributions, which have non-zero extent, the f 's have no simple dependence 
upon R, and except in the limits R -* 0 and R - ®. When I represents 

a nuclear attraction integral, only terms with is = 0 occur in the exjiansion. 
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In the form given by Eq. (2) the evaluation of an integral consists 

of determining the • Hence, the accuracy of an approximation for 

*1*2 N 

an integral may be evaluated in terms of the accuracy with which it approximates 
the • Therefore, predictions and evaluations of the accuracies 

I of phe Margenau approximations and the accuracies of the Sklar and MuUiken 
a^roximations, to which they are compared, will be carried out by considering; 

(1) the degree and manner in which the approximations include each f (and hence 
the degree to which the various multipole interactions are included) suad 

(2) the degree to which each f includes the effects of the overlap of the 
charge distributions. 


IV. THE MAR®NAU APPROXIMATIONS 

5 

The Margenau approximations vere employed on Intuitive grounds In cases 

vhere the ratios ^^^cd large* They were used to estimate 

and calculate values for the three* end four-center Integrals vhlch occurred 

In a determination of the H-Hg and H 2 -H 2 Intermole culetr potentleil energies at 

5-7 

Intermediate and long range distances. The results are In good agreement 

with calculations hy Mason apd Hlrschfelder for the same Intermoleculax 

10 

potential energies at short and Intermediate range distances. 

The Margenau approximations are apparently the original members of a class 
of approximations called the "point charis^e approximation."^ These 
approximations involve the replacement of the charge distributions In the 
Integreds by point charges. Different forms of the approximations Involve 
differences In the number of point charges per charge distribution^ differences 
In the positioning of the point charges, suid whether one or both charge 
distributions are so replaced. In all cases a distribution is replaced by 
point charges with the same total cheurge. Under the condition that a distribution 
is replaced by a single charge, the Margenau approximations are probably the 
most accurate of the point charge approximations that could be generally applied 
to the types of Integrals considered in Sec. IV> 

In the first and extremely simple approximation, Mj., Mtaurgenau re]^cea 
both cheo-ge distributions In Eq. (1) by point charges at the geometric centers 
of the dlatrlbutlona. Each of the charges has a value equal to the total 
charge of the distribution replaced. For the Integrals considered In Sec. VI 
the geometric centers are the points al> and In M^, therefore. 
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where the total charges Qx and 02 of "the two distrihutions of charge are 
equivalent to e times the overlap integrals, and A 2 , given hy 

A, » awd 

Equation (3) results from Eq . (2) v n only the Xx ® •*2 ® 0 term is 
taken and f is approximated hy its value in the limit ft/a^> H/'R., “• 

000 O GIO Cu 

Thus, Ml includes only the monopole-molopole term of an integreJ. end neglects 

the effects of overlap entirely. Mi wi;: > . xerefore he accuiate when 

ft/a^, and are simultaneously large but will be inewjcurate when 

R/a^ is small irrespective of the values of and ’^n R/a^ ■’ 0 

fiuid either and or R/Rq^^ and R/R^^ are fixed, the Mi approximations to 

the Integrals diverge, whereas the integrals themselves do not. 

In the limit R/a -• ® with R/R . and R/r , fixed. Mi does not in genereU. 
o ab cd 

yield a value for an integral that approaches the ■ xact value of that integral. 
In this limit eui integral as r-pproxlmated by Mj will be proportional to the 
exact result in the case in which the integral is one whose exact value 
approaches const, x QiQa/R, but an integral as approximated by hi will tend 
to zero relative to the exact result in the other case in vdiich the exact 
resTilt approaches const, x 4i%2 r/R* 

La Margenau’ s second approximation, Ma, the charge distribution of elec- 
tron 1 alcme is replaced by an equal point charge at the geometric center of 
*-he charge distribution, and the integration over the coordinates of election 
2 in Eq. (l) is then carried out. This process is repeated with the charge 
iistrlbution of electron 2 alone being similarly replaced, and the mean of 
the two results is taken. Thus, In M 2 . 
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I (a«,b^;cMO • 

♦ A,i«^Y„(?)VtIiUr,) , (5) 

where rg ' is the distance from electron 2 to the geometric center of the charge 
distribution of electron 1 and rj ' is analogously defined. Use of Mg therefore 
reduces the determination of a three- or four-center integral involving two 
charge distributions of non-zero extent to the determination of two three- 
center integrals that are equivalent to nuclear atti’action Integrals. For 
the cases considered in Sec. VI these three-center integrEJ.s are known in 
relatively simple, closed form. 

Equation (5) results from Eq. (2) in the following way. In Eq. (2) only 
■ 0 terms are taken and each f is given the form it would have if the 

OXgO 

charge distribution of flectron 1 were a point charge with a charge Qx. 

This process is repeated, taking only the • 0 terms, and the mean of the 

two results is teJeen. If the terms in the final result are then put in the 

form of Eq. (2), the resulting f is exact in the limit d/a , r/r , 

-* and the resulting for ^ i possess one-half 

of their exact values in this limit. Since only one charge distribution at 

a time is replaced by a point charge in Mg, each of the f' s includes the 

effects of overlap to a certain degree. Thus, Ma tedtes account of the 
mocopole-monopole term and a partial account of the monopble-hlgher multipole 
terms but neglects the higher multipole-higher multlpble terms, and for each 
of the included terms, Ms makes a partial allowance for ths effects of over- 
lap. Ms is therefore expected to be accurate for simultaneously large values 
of •'naccurate fox small R/a^, as v&a the case 
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vlth . Far aost Instances of tbe integrals, however, M2 vlU give more 
accurate values then because M2 Includes part of the effects of hl^ier 
nultlpdes and overlap idkereas does not. 

Li general, M2 like does not yield exact values for the integrals 
in tlK Halt ft/a^ « vlth ^ ^^ed 3°’*'*''*^* ^ aost, if not 

all eases in this lialt Mg gives axpressioas for the three- end fcur-cssuer 
integrals that are porapoortiaol to the eaaet expresaiOBS, ocnst. x ^Qg/M 
or ecnst* x Q^Qg in ft/g. 
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V. THE SKLAR AM) t«JLLIKEN APPROXD^IATIONS, 

' CCMPARISOMS TO Ml AM) Ma 

The accuracies of the Margenau approximations eure compared in Sec. VI 

8 0 

to the accuracies of the Sklar and MulUken*' approximations. The latter 

are chosen for this purpose, in spite of the fact that more accurate ap~ 

12 

proxlmations exist, hecause they yield relatively slmj^e formulae for the 
Integrals and hecause they have been the most extensively employed. 

the Sklar approximation, S, as applied to the cases of Sec. VI, the 
dbarge distributions In Eq. (l), e ® 

are replaced by charge distributions of equal total charge, e ^ * 

^j^p(ri) and e resulting charge distributions 

are centered at the points a& and Hence, S reduces three- and four- 
center integrals to two center integrals. In S, 

The Sklar approximation results from Eq. (2) vhen only the 0 

term Is retained and a particular approximate form for f^^^ Is employed. 

The approximate form of f^^ Is exact In the limit ^^^cd “* ** 

and Includes the effects of overlap. S Is therefore similar to in that 
It includes only the mocopole-monopole teirm of an Integral but Is superior 
to M]^, In particular for small ft/a^. In that It Includes the effects of 
overlap. Mg Is superior to S In that It takes partial account of the hl^^r 
multipole terms but Is inferior in that It takes less account of overlap. 

It is therefore expected that S will be more accurate than Ma for small R/a , 

c 

but Mg is expected to be more accurate than S for simultaneously large values 
of ft/a^, and r/R^^. 
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In the limit -» ® with fixed, values of integrals 

as approximated hy S approach the values of corresponding integrals as 
approximated hy . Hence, S will also yield values for the integrals 
that are sonetlmes proportional to the exact values in this limit. It is 
therefore eiqpected that the accuracy of S will usually he less than that of 
VLs hut greater than that of (since S includes the effect of overlap) for 
simultaneously large values of ft/a^, ^cd^®‘o* ^ limit 

ft/a^ -• 0 with ft/Rgjj charge distributions of S ap- 

proach the exact change dlstrlhutlons . Hence, the values of integrals 
given hy S are exact in this limit, and S is expected to he more accurate 
than Ml or M 2 for simultaneously small values of ft/a^, ^cd^*^o* 

In the MuUiken approximation, M, the charge dis trihut ions, e * (ri) 

' fllQf 

i|i^p(ri) and e irQ^(?a) iii Eq. (l) are replaced hy the charge dis- 

trihutiona of eqpal total charges, ^ 

§ ^cd^^cy^^^^* charge distrihution of electron 1 thus he- 

comes the sum of two charge distrihutions centered at nuclear centers a and 
h, and the charge distrihution of electron 2 hecomes the sum of two charge 
distributions centered at nuclear centers c and d. Hence, M reduCes tbree- 
and four-center integrals to a finite sum of two center integrals. In M, 
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The Mulllken approximation results frcsn Eq,. (2) vhen all terns in the 

multlpole expansion sure retained su^d peurticular approximate forms for the 

f. . are employed. Each of the approximate f^ . will Include the 
iiXsm sxism 

effects of overlap, hut each will in general he larger than the exact 

f^ . because M elongates the dis trihut ions along the intemuclear axes, 
iiiam 

Since M includes the effects of overlap to a larger degree than or Mg, 
it is e3(pected that and Hz will he less accurate than M for small values 
of RAq* Without a more detailed* consideration it is difficult to predict 
the accuracy of and Mj. relative to M for large values of 
and 

In the limit ft/a^ -* <» with aod fixed, values of the integrals 

as approriJaated l>jr M approach const. , ^ ^ ^ 

h are conincident with either of the centers c and d.due to an overemphasis 
of overlap at the coincident centers, or the values apf-rcach const, x 
otherwise. Hence, the accuracies of M may he inferior or superior to the 
accuracies of hut will usxially he inferior to the accuracies of Mg for 
simultaneously large values of ®ah^*^o* ^cd^®o’ true with 

S, the charge distrlhutlons of M approach the exact charge dletrlhutlons in 
the limit ■* 0 with fixed. Hence, M is expected to he 

more accurate than or M 2 for simultaneously small values of R/a^, 
and ^> 0 * 
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VI. DIRECT DETERMINATIONS OF AND Mg ACCURACIES, 

CCKPARISCNS WITH M AND S ACCURACIES 
The Integrals for which compariscaas between approximate and accvrate 
values were made were chosen from the limited number of relatively simple 
three- and four-center integrals for which accurate values are readily 
available. The angular configurations of the nuclear centers for each basic 
type of integral were taken to be as different as possible, and for each con- 
figuration all values of ft, R^^, and/or R^^ for which accurate values could 
be obtained were included, Hie formulae for the specific integreds in the 
nine subsequent cases eure given below. The nuclear configurations for each 
case are given in Fig. 2 . In case 5 'the nuclear centers lie in a single pl a n e , 
Cases 1 , 2, 5 , 4; 




rral 






dX 


( 8 ) 


Case 5‘ 





( 9 ) 


Cases 6, 7: 





W'WO* X 


e*''* 


i/*. 




( 10 ) 


Cases 8, 9: 






( 11 ) 
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In cases 1-5 the sequences of configurations are of the first type in 

which ft/a^ varies while and are fixed. In cases 6-9 the sequences 

are of the second type in which ft/a^ varies while *^^^cd fixed. 

Mg does not apply to the three-center nuclear attraction integrals of cases 

l8 

1’4 whose evaluated form is known because one of the two terms that are 
averaged in Mg Is the exact result in this case. Cases 6-9* aod to a lesser 
degree the other cases, represent e severe test for Mx and Mg since these 
cases largely involve relatively small values of and/or whereas 

the Intent of and Ms is to he accurate tor large values of R/h^'^ 

R/I^cd* 

The Mx, .is, S, and M foraiulae for the Integrals of each case may he 
obtained from the equations and discussions of Sec. II- IV. For M^, S, and 
M determination of the integrated formulae involve oiily simple, well known 
two-center integrations. Ih Ms, a three-center nucleeu: attraction Integral 
occurs, the evaluated form of which Is given in Ref. l8. 

The results of the comparisons are given in Tables 1-9 for cases 1 - 9 * 

The accurate values of the integrals are tedcen from the works of Hirschfelder 

19 3 k 

and Weygant, Boys and Shavltt,^ and Magnasco and Musso. 

In considering the results it is important to bear in mind that neither 

Ml, Ms, S, nor M yield values ,for the integrals that are in general either 

upper or lower bounds to the exact values. Thus, t>e difference between 

accurate and approximate val'^ each c«4e may pass through zero within 

the particular sequence of configurations of each case. In lieu of means by 

vrhich these points of zero error can be predicted, the high accuracy in the 

vicinity of these points must be regarded as fortuitous and should not be 

considered to have any particular significance in assessing the accuracy of the 

approximations . 
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The exact R = ® limiting forms could not be determined for cases 
8 and 9* The R « <» "accurate" ved.ue for case 8 was found by replacing the 
charge distributions by uniform line segments of charge extending from a to 
b jtnd c to d. This process gives a result which is proportional to the exact 
value but gives only an approximate value for the constant of proportionality. 
For case 9 the above model fails and only an expression which is proportional 
to the exact result could be determined. 

As expected, it is /-.pparent iVom the tables that Mx possesses poor ac- 
curacy in the majority of instances considered here and is usually inferior 
to Mg, S, and M. In most cases, therefore, is useful only in obtaining 

very rough estimates for the values of the Integrals. Exceptions to this 

occur, as expected, for simultaneously large values of ft/a , r/r . , and ft/R ^ 

O ED ' cd 

where is nearly as accurate as Mg axxl B and more accurate than M, as il- 
lustrated by Tables 1-^. Tables 1-^ and to a lesser degree Tables 6-9 
demonstrate the fact that as ft -» <» values given by and S become equal. 

Since S in cl u des the effects of overl«ip to a large degree while negpLects 

overlap entirely, values of ft for which Mi and S give nearly equal results 
marF the points at which the effects of overlap become unimportant. Such 
values of ft lie between 5 suid 7 Tor Cases 1-5 and are estimated to be 
about the same for Cases 6-8. The value for Case 9 is uncertain. Due to 
the simplicity of Mi, It should be useful in determining which, if any, 
integrals may be ne^ected in a given molecular structure calculation. 

The results given in all of the tables bear out to a large degree the 
expectations regarding the accuracy of Mg relative to S and M expressed in 
Sec. V. A comparison of the accuracies of Mg, S, and M for sequences of 
values of ft in which R^^ and R^^ remain fixed is illustrated only by Table 5. 
Here Mg is less accurate than S and M for small values of ft but more accurate 
than S for ft 5 and more accxirate than M for A 3 * 5*5 Compeurisons of 
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accuracies for sequ ences In vhich remain fixed occvir in Tables 

6 - 9 . Here the maximum values of fi included in the tables (except for fi -» •) 

vary frcxa about a for Table 6 down to about I .5 for Table 9* It is 

therefore not unexpected that in these cases % is inferior to both S and M 

in nearly all instances. It is clear, however, from the values and the 

tendencies of the values of the last rows of the tables and from the rt -* oe> 

limiting values that M 2 becomes more accxirate than S and M for certain values 

of R, as was expected. For Table 7 these values of R are about 4 and 3 a^ 

relative to S and M respectively, and for Table 8 the value relative to M is 

about 3 a. ]n the other cases these values of R lie off the tables but 
o 

appear to be roughly the same. It may be concluded, therefore, that in 

nearly all cases M 2 is more accurate than S for values of R 5 a^ and more 

acc\irate than M for R 3*5 a^. An exception to this would probably occur 

for cases not included here in which r/R . and r/R . are considerably smaller 

ab cu. 

than unity. 

It is Interesting to note that in all cases considered here S is more 

accxirate than M 2 for R < 4 a and that S and M are more accurate than one 

~ o 

another for the same values of R in about as many Instances as not* Since S 
includes the effects of overlap but does not Include the effects of higher 
multi poles, while Mg and M Include both effects. It may be concluded that 
the effects of overlap are more Important than the effects of hlg^r multlpoles 
In most cases for A < ^ a^* Since Mg Is more accurate than S for A 9 
the opposite Is true for these higher val\ies of R. This is consistenc with 
the conclusion previously reached by comparing and S that the effects of 
overlap are mimportant for R > 7 a . 

•N# O 
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20 

Ihe statement has been made that the smallest of two approximate 
valiies for an Integred. Is usually the more accurate. comparing the 
Sklar and 14 ulliken approximations for the cases considered here it is 
found that there are 36 Instances in which the smaller value is more 
accurate and I9 Instances in which it is not. 
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VII. A CCMBUtATION OF THE MARGENAU APPROXIMATIONS 


An approximation that is often superior to both and Mg may be obtained 
by taking twice the value of an integral as gi\'en by Mg and subtracting from 
that the value as given by Mj . There are two reasons why this should be an 
improved approximation. First, it may be seen from the discussion of Sec. IV 
that in the limit R -» » with R^^ and R^^ fixed, 2Mg-Mi gives exactly the monopole 
monopole and monopole -hi^er multipole terms of the integral. This is in 
contrast to Mj which includes only the monopole -monopole term and to Mg which 
in addition includes only cne half of each monopolerhigher multipole term. 
Second, 2Mg-M^ will represent the effects of overlap to a better degree for 
values of R that are not very small. This is because the partial allowance 
of overlap in Mg is improved upon by the subtraction of Mj which includes no 
effects of overlap. In other words, the lack of complete allov^nce for the 
effects of overlap in Mg is in some degree canceled by subtracting a quantity 
which has a greater lack of allowance for overlap. For very small values of 
R, 2Mg-Mi will be poor due to the divergence of M^ for R 0. 

Values of the Integrals approximated by ^g-M^ fo. Cases ^-9 are given 
in Tables ^- 9 . Case ^ is the only case in this group which involves a sequence 
of configurations In which R varies while and remain fixed. Txx Case ^ 
accurate than and Mg for R > a^ is more accurate than S and 
M for R > }a^, is at least neaxly as accurate as S and M for R a^, and is 
less accurate than Mg, 3, and M for fi <»■• 

For Cases 6 and 7, where H^^/r rema±3s fixed as R is varied, 2Mg-M^ is 
exact in the limit R -» •. This is due to the :teet that one of the charge 
distrlbutlcns is spherically symmetric, and the integrals therefore do not 
contain any of the higher multlpole-higher multipole terms which are neglected 
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by 2Ms“Mi. In spite of the fact that increases in proportion to R as 

ft ”♦ », 2 H 2 “Mi gives exactly each monopole-multipole term (each approaching 

eA X ft”^ as ft “* “) in this limit. In Case 6, 2 Ml”Mi is superior to the other 

approximations for ft > 2a^. In Case T, where the effects of overlap are 

greater due to the different nuclear configuration, is superior to M 

for ft > superior to S for ft > 4a^, and comparable or superior to 

for ft > 2a . 

o 

In Case 8, where and are fixed as ft varies, the integral 

includes higher multipole-higher multipole terms emd hence 2Ife-Mi is not 
exact in the limit ft -» ®. The limiting value is nevertheless better than 
that given by the other approximations. Here 2M2-Mi is superior to M for 
ft > 3a , superior to S for > a vedue of (probably) about 4a , and ccanparable 
or superior to Ma for ft >^a^. In Case 9> where and are fixed 

6uad where lines ab and bd partially overlap, 2 M 2 -Mi like Mg has the proper 
form is the limit ft “* ® whereas the other approximations do not. Here 
2Ma-Mi is apparently rou^ly cca5>arable or superior to Mg, S, and M for 
ft ^ 1*^ a_. 

It may be concluded, therefore, that in most cases 21fe-Mi Is more 
accurate than S for ft > 4a , is more accurate than M for ft > 3a , emd has 

ew O O 

about the same or better accuracy than Ife for ft > 2a^. Exceptions to these 
conclusions might be coses involving vedues of ft con^irable to those considered 
here but involving much larger values of Rjj^ than consldeied. 
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VIII. SUMMARY AND COMMENTS 


The Margenau app_oximations to three- and four-center moleculeu: integrals 

have heen tested for accuracy and compared with the Sklar and MuUiken ap- 

-r/a 

proximations for many cases involving atomic orbitals of the form e ’ o. The 
former approximaticos were designed to be accurate for cases involving large 
distances between charge centers . It was found that Margehau* s first very 
simple approximation is accurate for cases in which the distance between the 
centers of charge is at least about five times as large as a^ and at the same 
time the intemuclear distances of each charge distribution ears less than or 
about equal to a^. In such cases this approximation is nearly as accurate 
as the Shlar approximation and is more accurate than the MuUiken approxima- 
tion. Meurgenau' s second, more complicated approximation was found to be 
accurate in the same cases and for these cases is more accurate than the other 
approximations considered. !Ihls approximation was also found to be more 
accvirate than the other approximations vhen the intercharge separation is 
again at least about ^ a^ but where the intemuclear distances have any value 
except values that are very large relative to the intercheurge distance. 

An approximation was proposed in which an integral is taken equal to 
twice the value given by Margenau' s second approximation minus the value 
given by Margenau' s first approximation. This approximation was found to 
have an accuracy about equal to or better than Meirgenau's second approxima- 
tion in all instances except those in which the Intercharge sepeiratlon is 
small. Furthermore, the approximation was found to be superior to the Sklar 
and MuUiken appio:clmations fax* somewhat smaUer values of intercharge 
separation than was true of Margenau' s second approximation. 

The results obtained here may be generalized to cases where the basic 
atonic orbital is e ' ° by an appropriate change of distance scale. 
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In addition^ a means vas presented for predicting and \ 2 nderstandlng the 
accuracies of molecular Integral approximations. In vhlch consideration Is 
given to tbt effects of overlap and higher nultlpoles at the charge distri- 
butions. It vas found that predictions of the relative accuracies of the 
approximations considered here vere borne out. This means of considering 
the apj^oximatlons led to the combined form of the Margemu approximations. 

For small values of ft/a^ the accuracy of any of the approximations con- 
sidered here is prctbably Izfisufflclent for all molecular structure calcula- 
tions, althou^ they may be adequate for determining the conflguratlaa of 
the molecule. Approximations ^dilch are sometimes more accurate than those 
considered here have been proposed and eoqlciyed. There have been many sug- 
gested generalizations of the Sklar Kulllken approximations vhlch 

17 12 

Improve their accuracy to some degree. The Clzek approximations oifer 

a more substantial lBq;irovmBent. In the latter each chaz^ distribution Is 

replaced by a sum of sj^ierlcally sympietric charge distributions placed aXong 

the intemudeeur eixls in such a fashion that the first few iz terms of 

Eq. (2) are exact In the limit Jl/a^ -• ® with and fixed. In this 

process the effects of overlap are apparently veil accounted for. The 

accuracies for conflguratloa sequences in which R varies idille ^^^cd 

ronaln fixed (the only case considered by Clzek) are quite good, and the 

* 

accuraciea for a I m u ltaneously large values of ft/a^, ^ 

good. However, the approximations would be expected to be less accurate for 

simultaneously small values of ^ latter case 

improved accuracy might result if the placement of the approximate charge 

distributions were such as to yield the first few fi, jlz terms of Eq. (2) 

exactly in the limit ft/a -• 0 with R . /a and R ,/a^ fixed, 

' ' o ab o cd' o 
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17 

Hirschi'eHer has made direct use of Eq. (2) with approximations to 

1 IP 

the first few ^ ' Silverstone and Rudenberg * have used Eq. (2) 
but with the bipolar expansion centers placed on the nuclei, to obtain 
exact, higtj-y complicated expressions for some of the three- and four- 
centar Integrals. This author has given consideration to obtaining exact 
expressioocs for each of the ^ ®q. (2) for the case in which the 

expB:'uion points are chosen to be the geometric charge centers ab and cd 
in <-.rd ir to obtain more rapid convergence of the e^gpansion. At the present 
time an exact multipole expansion for the charge distribution eCna^)"*^ 
about the point ab has been obtained. 
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FIGURE CAPTIONS 


Fig. 1. Hie varlovis line segments, distances, angles and particles involved 
in the discussions of the text. The nuclei are particles a, b, c, 
and d and the electrons are particles 1 and 2. Points db and cd 
bisect line segments ab and cd respectively, and particles and points 
a, b, c3., P, and a' lie in a single plane, as do c, d, ab, P, emd c' . 
The self-descriptive distances r^^, r^^^, r^^ and r^^^ are not shown. 

Fig. 2. The nuclear configurations involved in the integrals of cases 1-9 . 

All nuclei may be considered to lie in the plane of this page. Dots 
mark the midpoints of the line segments on idiich they are placed. 
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p 

1 

1 

accurate 

(ref.) 


S 

M 

1 

0 

0.782455 

(19) 

OO ' 

(+ 00%) 

j 

0.858385 

(+9.70%) 

0.769424 

(-1.66%) 

3 

1 

- i 

1.4142 

0.509687 

a9) 

0.606970 

(<^19.1%) 

0.520359 

1 

(+2.09%) 

0.501029 

(-1.70%) 

5 

2.4495 

0.339261 

(19) 

0.350434 

(+3.29%) 

0.341424 j 

(+0.638%) 

0.335257 

(-1.18%) 


3.4641 

0.246089 

(19) 

0.247794 

(+0.693%) 

0.246711 

(+0.253%) 

0.244247 

(-0.743%) 

9 

4.4721 

0.191549 

(19) 

0.191941 

(+0.205%) 

0.191804 

(+0.133%) 

0.190623 

(-0.483%) 

13 

6.4807 

0.132370 

(19) 

j 

i 

0.132452 

(+0.062%) 

0.132450 

1 (+0.060%) 

• 

0.132057 

(-0.236%) 

no 

e« 

0.858385/R' 

0.858385/%' 

(0%) 

0.858385/%' 

(0%) 

0.858385/%' 

(0%) 


Table 1. Comparlaon of aoourato aitd i^tpgoarimate valuea for tbe tntesral 
of case 1. Values of the inte^pral have been divided by e > R ■ a^ . 

P = 1^1 + (2^iyR)*. SP = S/a„. 
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p 


accurate 

Orel.) 

Mi 

8 

M 

1 

0 

0.464605 

ao) 

oo 

<^%) 

0.586453 

(+26.2%) 

0.427717 

(-7.94%) 

3 

2.8284 

0.201223 

(19) 

0.207342 

(h3.049^ 

0.204569 

(+1.66%) 

0.193546 

(-3.82%) 

5 

4.8990 

0.119160 

a&> 

0.118709 

(+0.461%) 

0.119670 

(H>.430%) 

0.117259 

(-1.60%) 

7 

6.9282 

0.084471 

a9) 

0.084647 

(H).208^ 

0.084647 

(+0.208%) 

0.083778 

(-.820%) 

9 

8.9443 

0.065486 

ao) 

0.065567 

<+0.124%) 

0.065567 

(fO.124%) 

0.065161 

(-0.496%) 

m 

•• 

0.858385/Sr 

0.858385/S(’ 

(0%) 

0.858385/%' 

(0%) 

0.858385/%' 

(0%) 


Table 2 . Con^iariSKMi of accurate and iq>proirimate vab&ea for the integral 
of case 2. Values of the Integral have been divided by e^a^. R = 2a^. 

P * 1^1 + (2*/R)a. *• » */».. 
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p 

X/a„ 

accurate 

(ref.) 

Mj 

S 

M 

0 

0 

0.782455 

(19) 

00 

(+ 00 %) 

0.858385 

(+9.70%) 

0.7o9424 

(-1.66%) 

1 

H 

0.735758 

(19) 

1.716771 

(+133.3%) 

0.769424 

(+4.58%) 

0.742216 

(+0.878%) 

3 

1.5 

0.500574 

(19) 

0.572257 

(*■14.3%) 

0.501029 

(+0.091%) 

0.515828 

(+1.05%) 

5 

2.5 

0.336648 

(19) 

0.343354 

(+1.99%) 

0.335257 

(-0.413%) 

0.344451 

(+2.32%) 

7 

3.5 

0.245086 

(19) 

0.245253 

(+0.068%) 

0.244247 

(-0.342%) 

0.248764 

(+1.50%) 

9 

4.5 

0.191074 

(19) 

0.190752 

(-0.168%) 

0.190623 

(-0.236%) 

0.192933 

(+0.973%) 

13 

6.5 

0.132212 

(19) 

0.132059 

(-0.116%) 

0.132057 

(-0.117%) 

0.132842 

(+0.476%) 

00 

00 

0.858385/X' 

0.858385/9r 

(9X) 

0.858385/31' 

(0%) 

0.858385/9C 

(0%) 


Tabi« 3. CoxDiMurison of accurate and approximate values for the Integral 
of case 3. Values of the Integral have been divided by eVa^. R a^. 
p * 2«/a^. r » 
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p 

*/a. 

accui'ate 

(ref.) 

M, 

S 

M 

0 

0 

0.464605 

(19) 

OO 

(+»»%) 

0.586433 

(+26.2%) 

0.427717 

(-7.94%) 

1 

1 

i 

1 i 

0.406005 

(19) 

0.586453 

(+44.4%) 

0.427717 
(+5.35%) \ 
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(+6.35%) 

3 

3 
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0.193546 
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0.211741 

(+7.66%) 

5 

5 
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(+3.27%) 

7 

7 
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(-0.417%) 

0.083778 

(-0.418%) 

0.085522 

(+1.65%) 

9 

I 

9 

0.065324 

(19) 

0.065161 

(-0.250%) 

0.065161 

(-0.250%) 

0.066278 

(+0.998%) 

ao 

oo 

0.858385/S’ 

0.858385/3C 

(0%) 

0.858385/S’ 

(0%) 

0.8S8385/S(‘ 

(0%) 


Table 4. Comparison of accurate and approximate values for the int^^al 
of case 4. Values of the integral have bemi divided by eVa^. R = 2a^. 

P = 2»/R. = »/a^. 
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Table 5. CcaiQMrlSQA ot accurate and approximate ralues for the integral of case 5. Values of the integral 
have been divided by eVa^. R = l.t»9a^. 3t‘ = 3t/a^. 


accurate 
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have been divided by eVa^. S’ - |S R*/2. A = (R*V3 + R* + 1) eip(-R'), 
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Table 7. Comparison of accurate and approximate values for the integral of oaae 7. Values of the integral 
have been divided by eya,. JR* = 3R’/2. A « ^S’Vs •♦• H' + l)exp(-R’),R = R/a . 
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Table 8. Comparison of accurate and approximate values for the Integral of case 8. Values of the integral 
have been divided by eVa«. 31 = R . A * (R'Vs + R + 1) exp(-R'), R' = R/a^. 31' = 31 /a . 
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Table 9. Comparison of accurate and approximate values for the integral of case 9. Values of the integral 
have been divided by e^/a,. R /2. A. » (H'Vs + R' + 1) exp(-R'), A. » (4R'V3 + 2R' + 1) exp(<2R'). R* = ^ 
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